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ABSTRACT 

Context. Thanks to its excellent 5100 m high site in Chajnantor, the Atacama Pathfinder Experiment (APEX) systematically explores 
the southern sky at submillimeter wavelengths, in both continuum and spectral line emission. Studying continuum emission from 
interstellar dust is essential to locating the highest density regions in the interstellar medium, and deriving their masses, column 
densities, density structures, and large-scale morphologies. In particular, the early stages of (massive) star formation remain poorly 
understood, mainly because only small samples of high-mass proto-stellar or young stellar objects have been studied in detail so far. 
Aims. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand 
how and under what conditions star formation takes place. Only a systematic survey of the Galactic Plane can provide the statistical 
basis for unbiased studies. A well characterized sample of Galactic star-forming sites will deliver an evolutionary sequence and a 
mass function of high-mass, star-forming clumps. This systematic survey at submillimeter wavelengths also represents a preparatory 
work for Herschel and ALMA. 

Methods. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LAB OCA) 
is a 295-element bolometer array observing at 870 yt/m, with a beam size of 19'.'2. Taking advantage of its large field of view (11'4) 
and excellent sensitivity, we started an unbiased survey of the entire Galactic Plane accessible to APEX, with a typical noise level of 
50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). 

Results. As a first step, we covered ~95 deg^ of the Galactic Plane. These data reveal ~6000 compact sources brighter than 0.25 Jy, or 
63 sources per square degree, as well as extended structures, many of them filamentary. About two thirds of the compact sources have 
no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. 
Other compact sources harbor hot cores, compact Hii regions, or young embedded clusters, thus tracing more evolved stages after 
massive stars have formed. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 
10 to a few 100 M©. In this first introductory paper, we show preliminary results from these ongoing observations, and discuss the 
mid- and long-term perspectives of the survey. 
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1. Introduction 

Dust continuum emission in the (sub)millimeter (submm) range 
is one of the most reliable tracers of the earliest phases of star 
formation since it directly probes the dense interstellar material 
from which the stars form. Trying to understand the formation 
and early evolution of stars is an important field of mo dern as- 
trophysics (see the review o f iMcKee & Qstrikedl2007L and ref- 
erences therein). 
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A large amount of theoretical effort has been aimed at con- 
straining the early stages of the formation of isolated low mass 
(solar-like) stars, although two quite different pictures are be- 
ing maintained. On the one hand, a collapse picture has emerged 
whose "long" timescales are determ ined by ambipolar diffusion 
(IShu et al.lll987l:lAdams et al.lll987h . On the other, a scenario is 
advocated in which the star formation rat e is governed by su- 
personic turbulence (iMac Low & Klessenll2QQ4l) . Many obser- 
vational campaigns have contributed crucial data and in partic- 
ular systematic (sub)millimeter surveys of nearby star-forming 
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regions, such as Chamaeleon, Perseus, and Ophiuchus, have 
resulted in the detection of many prestellar and proto- stellar 
condensations ( Andre et al. 2000; Ward-Thompson et al. 2007). 
They have allowed determination of the dense core mass spec- 
trum down to substellar masses, and investigation of its relation- 
ship to the initial mass fun ction (iMotte et al.ll 19981: iGahm et aP 
l2QQ2l:lHatchell et al.l l2005^. 

High-mass star forma t ion is even less well constrained (e.g., 
IZinnecker & Yorkd[2QQ7l: iBeuther et al.ll2007L for reviews) and 
there is controversy about how these stars at tain their fi nal 
masses: either driven by tu rbulence jM cKee & Tad l2003l) or 
competitive accretion (e.g., iBonnell et al.l ^2004). Since high- 
mass stars seem to form mostly in clusters (along with many 
more lower-mass stars), it is fundamental to have a good esti- 
mate of the mass distribution of cluster-forming clumps in the 
Giant Molecular Clouds (GMCs) they are part of. This is exactly 
one of the major goals of the present eff'ort. 

A cardinal reason for the many open questions concerning 
high-mass star forma tion is that the relevant timescales are short 
(iMcKee & Tanll2002l) , implying low number statistics in finding 
stars in a given evolutionary state. Moreover, high-mass stars 
are rare. This implies that regions in which they form are gen- 
erally at large distances, which leads to the basic problem that 
there objects are difficult to identify. For instance, little is known 
about the Galactic molecular ring, a few kpc distant from us, 
where most of the st ar formation is presently taking place (e.g., 
iBronfman et aLlf200Q) . 

While various samples of high-mass proto-stellar objects 
have been defin ed based on various criteria, such as maser 
emission (e.g., Walsh et al. 1997'; "Pestalozzi et al.' '2005) and 
infrared (IR) colors (Pallaetal. 1991: Molin ari et al. 1 996 



ISridharanet aP 120021: iLumsdenet al 
120081), an unbiased sample has not yet been assembled 
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effort has focused on follow-up studies of bright IR sources de- 
tected by the Infrared Astronomy Satellite (IRAS) and show- 
ing fa r-infrared (FIR) colors typ ical of ultra-compact Hii re- 
gions (IWood & Churchwellll 19891). e.g., in high-density m olec- 
ular tracers feronfman et al." "1996; "MolinanetaLl ll996|), dust 
continuum emission (Faundez et al. 2004), or maser emission 
(e.g.. iPalla et al.lll993l: IBeuther et al.ll2002D . While some stud- 
ies have revealed a few massive cold c ores in the neighbo rhood 
of the central compact Hii region re.g. JGaravet ZI I2004I) . they 
were mostly biased against the earliest, possibly coldest phases 
of massive star formation. 

Alternati vely, dense conden sations in Infrared Dark Clouds 
(IRDCs, see ISimon et al.|[2006i. for a compilation) are promis- 
ing hunting gr ounds for the initial s tages of high-mass star 
formation (e.g ., iRathborne et al. 2006; see also the reviews by 
iMenten et all 1200 5: Wyrowski 2008). IRDCs appear in absorp- 
tion by dust even at mid-infrared wavelengths against the dif- 
fuse emission of the Galactic Plane and trace high column den- 
sities. However, IRDCs become increasingly difficult to identify 
the further away they are, and the relative contributions by fore- 
ground and background emission may be difficult to estimate. 
Moreover, not all massive dust condensations appear as infrared 
dark clouds. For example, the two most massive condensations 
(with masses of 200 and 800 M© at a distance of 1.3 kpc) ob- 
served at 870 yum in the photo-dissociation region (P DR) bor- 
dering the RCW 120 Hii region (D eharveng et al.ll200 9) are not 
detected as IRDCs at 8 or 24 jim. They are hidden, at these wave- 
lengths, by the bright emission of the adjacent foreground PDR. 

Cold dust absorbing the IR radiation also shows thermal 
emission in the submm regime. This grey body emission is gen- 
erally optically thin and, thus, an excellent tracer of the dust mass 



and, by inference, total mass of the emitting cloud. Therefore, 
observations in the submm continuum are complementary to the 
IR extinction techniques. Given the short timescales for the for- 
mation of high-mass stars, all early formation stages will still 
be associated with strong dust continuum emission, allowing a 
systematic study of a large range of evolutionary phases. 

Several groups have performed large-scale mapping of 
molecular complexes in the (sub)millimeter continuum. Using 
the Submill imetre Common User B olometer Array (SCUBA) 
instrument, iJohnstone et al.l (|2004|) mapped 4 deg^ o f the 
Ophiuchus star-forming cloud, and iHatchell et al.l (l2005l) cov- 
ered 3 deg^ of the Perseus molecular cloud. Complementary 
to th e Spitzer ''Cores to Disks" legacy project (Ev ans et al] 
l2003l) jEnoch et al.l (l20Q6l) m apped 7.5 deg^ covering the Perseus 
cloud, and I Young et al.l (12006) mapped nearly 11 deg^ covering 
Ophiuchus, both using the BOLOCAM instrument at 1.1 mm. 

In the field of high-mass star formation jMoore et al.l (l2007h 
covered 0.9 deg^ of the W3 giant mo l ecular cloud with SCUBA. 
Going to larger scales, iMotte et al.1 (l2007h used the MAMBO 
instrument to obtain a complete dust-continuum map of the 
molecular complex associated with the Cygnus X region. From 
this 3 deg^ map, they could derive some statistically signifi- 
cant results about the characteristics (e.g., mass, density, out- 
flow power, mid-IR flux) and lifetimes of different evolutionary 
phases of high-mass star formation. In particular, no high-mass 
analog of pre- stellar dense cores was found by this study, hint- 
ing at a short lifetime for the earliest phases of dense cores in 
which high-mass stars are being formed. An unbiased survey of 
the complete Galactic Plane will allow us to place far stronger 
constraints on this stage of high-mass star formation. 

On an excellent 5100 m high site and at a latitude of - 23°, 
the 12 m diameter APEX telescope (iGiisten et al.1 l2006allbh is 
ideally located to make sensitive observations of the two inner 
quadrants of the Galactic Plane. A consortium led by the Max 
Planck Institute fiir Radioastronomie (MPIfR) in Bonn, involv- 
ing scientists from the Max Planck Institute fiir Astronomic in 
Heidelberg and from the ESO and Chilean communities has em- 
barked on the APEX Telescope Large Area Survey of the Galaxy 
(ATLASGAL). This project aims at a systematic survey of the 
inner Galactic Plane, mapping several 100 deg^ with a uniform 
sensitivity. In this paper, we present the first data taken and an- 
alyzed. The observing strategy and data reduction are described 
in Sect. [2l and Sect. [3l First preliminary results are discussed in 
Sect, m Finally, we discuss longer term perspectives in Sect. \5\ 



2. Observations 

Observations of continuum emission from interstellar dust 
at (sub)millimeter wavelengths are completed most effec- 
tively with sensitive broadband bolometer detectors. Modern 
(sub)millimeter bolometer development started in the mid 1980s 
with single-pixel-element instruments, with a resolution, i.e., 
diffraction-limited FWHM beam size, of IV for the IRAM 
30 m telescope at 1.2 mm and 19'' for the James Clerk Maxwell 
Telescope (JCMT) at 1 . 1 mm. Since then, arrays with more and 
more elements have been developed and deployed to telescopes. 
Modern arrays of bolometers have hundreds of elements: 384 for 
the Nyquist- sampled Submillimeter Hig h-Angular Resolution 
Camera 2 (SHARC II. lDoweU et al.l2003h . and 295 for the Large 
APEX Bolometer Camera (LAB OCA), which was recently de- 
ployed at the APEX telescope ( Siringo et al. 2007, 2009). With a 
field of view of 1 1 f 4, and a single-pixel sensitivity (noise equiva- 
lent flux density, NEFD) in the range 40-70 mJy s^/^ LABOCA 
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Fig. 1. Coverage of the ATLASGAL observations overlaid on an IRAS false color image (12 fim in blue, 60 jim in green and 
100 yum in red). This image covers ±90° in galactic longitude, and ±10° in latitude. The long frame shows the area that we plan to 
cover until the end of 2009: ±60° in / and ±1.5° in b. The smaller frames inside delineate the area that was observed in 2007 (see 
also Fig. [3]). 



is the most powerful bolometer array for large-scale mapping 
operational on a ground-based telescope. 

The ATLASGAL observations are carried out at the APEX 
12 m submm antenna, which is a modified copy of the VERTEX 
prototype antenna for the Atacama Large Millimeter Array 
(ALMA). The antenna has a surface accuracy of 15 jjm rms, i.e., 
higher than ALMA specifications. Two additional Nasmyth cab- 
ins allow several receivers to be ready for use. The LAB OCA 
instrument is an array of 295 bolometers arranged in an 
hexagonal pattern, w ith two-beam spacing between bolometers 
dSiringo et al. 2009). It is located at the Cassegrain focus, where 
it has an eff'ective field of view of 1 1 f 4 in diameter. Small vari- 
ations in the beam shape are seen within the fiel d of view, but 
should not aff'ect the data (see lSiringo et al.ll20^9 . for a discus- 
sion of these eff'ets). Its bandpass is centered on 870 jim (fre- 
quency 345 GHz) with a bandwidth of 60 GHz. The beam at 
this wavelength was measured to have a width of 19^.^2 FWHM. 
The number of usable bolometers at the end of the reduction, 
after flagging the noisy ones and the ones that do not respond 
(Sect. [3]) is usually between 250 and 260. 

In the present paper, we focus on the data acquired in 2007, 
which cover 95 deg^: -30° < / < +11.5° and +15° < / < 
+21°, with \b\ < r (Fig. [B, at a one-cr sensitivity of 50 
to 70 mJy/beam. For comparison, the 9-year lifetime of the 
SCUBA instrument resulted in a total of 29.3 deg^ mapped at 
850 Jim (iDi Francesco et al.ll2008h . A single observation consists 
in an on-the-fly map, 2° long and 1° wide, with 1 f5 steps between 
individual lines, resulting in fully sampled maps covering 2 deg^. 
The scanning speed used is 37s. This fast scanning allows us to 
perform observations in total power mode, i.e., without chop- 
ping. Each map crosses the Galactic Plane with a ±15° position 
angle with respect to the Galactic latitude axis. Two consecutive 
maps are spaced by 0.5° along the Galactic Plane, with a 30° 
difl'erence in position angle (Fig.©. Thus, each position on the 
sky is observed in two independent maps, of difl'erent position 
angles, greatly reducing many systematic efl'ects (e.g., striping 
in the scanning direction). 

In the longer term, we plan to extend the coverage to ±60° 
in longitude and ±1.5° in latitude (Fig.O, with the same sensi- 
tivity of 50 mJy/beam. An ESO large programme was approved 
to conduct these observations, which will require a total of 430 
hours, including all overheads, to cover 360 deg^. The observ- 
ing time will be shared between 45 hours from Chile, 225 hours 
from the Max-Planck-Society (MPG), and 160 hours from ESO, 
spanning four semesters in 2008-2009. 



3. Data reduction and source extraction 

The raw data are recorded in MB -FITS (Multi-Beam FITS) for- 
mat by the FITS writer, which is pa rt of the APEX Control 
System (APECS, iMuders et al] l2006h . These data were pro- 
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Fig. 2. Mapping strategy: each rectangular frame represents a 
single map. The plus symbols show the maps centers. A position 
angle of 15° with respect to the Galactic axis was chosen. Two 
neighboring maps are observed with position angles of -hi 5° and 
-15°, respectively. 

cessed using the BOlometer array data Analysis package (BOA; 
Schuller et al. in prep.), which was specifically developed as part 
of the LAB OCA project to reduce data obtained with bolometers 
at the APEX telescope. The reduction steps involved in the pro- 
cessing of the ATLASGAL data are described in the next sub- 
sections. 

3.1. First step: from raw data to 2 deg^ maps 

Each single observation is processed and calibrated to compile 
a map covering 2 deg^. The data are calibrated by applying an 
opacity correction, as deter mined from skydips observed typ- 
ically every two hours (see ISiringo et al.l l2009l for more de- 
tails). In addition, the flux calibration is regularly (every hour) 
checked against primary calibrators (the planets Mars, Neptune 
and Uranus) or secondary calibrators (bright Galactic sources, 
for which the fluxes at 870 jim were measured during the com- 
missioning of LAB OCA). If any discrepancy is found between 
the opacity-corrected flux of a calibrator and its expected flux, 
then the correction factor required to obtain the expected flux is 
also applied to the science data. Given the uncertainties in the 
fluxes of the planets themselves, we estimate that, altogether, 
the typical calibration error should be lower than 15%. Pointing 
scans are also completed using bright sources every hour; the 
pointing rms is typically of order 4'\ 
The reduction steps are as follows: 

- flagging bad pixels (those that do not respond, or do not see 
the sky) 

- correlated noise removal (see below for details) 

- flagging noisy pixels 

- despiking 

- low-frequency filtering (see below) 
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- first-order baseline on the entire timestream 

- building a map, using natural weighting: each data point has 
a weight 1/rms^, where rms is the standard deviation of each 
pixel 

The main source of noise in the submm is variable atmo- 
spheric emission (sky noise), which varies slowly with typical 
amplitudes of several 100 Jy. Fortunately, this sky noise is highly 
correlated between bolometers. For each integration, the corre- 
lated noise is computed as the median value of all normalized 
signals. It is removed in successive steps: as a first step, the me- 
dian for all bolometers is computed, and then the median for all 
pixels connected to the same amplifier box (groups of 60 to 80 
pixels); finally, the correlated signal is computed for groups of 
20 to 25 pixels sharing the same read-out cable. A consequence 
of this processing is that uniformly extended emission is filtered 
out, since there is no way of distinguishing extended astronom- 
ical signal observable by all bolometers from correlated (e.g., 
atmospheric) noise. The strongest eff'ect occurs when computing 
the correlated noise in groups of pixels: these groups typically 
cover ~ 1^5 X 5' on the array, and since the correlated noise is 
computed as the median value of the signals, any emission that 
is seen by more than half the bolometers in a given group is sub- 
tracted. Thus, uniform emission on scales larger than ~ 2f5 is 
filtered out. 

Another low-frequency filtering is performed, to reduce the 
eff'ects of the slow instrumental drifts, which are uncorrected be- 
tween bolometers. This filtering is performed in the Fourier do- 
main, and applies to frequencies below 0.2 Hz, which, given the 
scanning speed of 37s, translates to spatial scales of above 15', 
i.e., larger than the LABOCA field of view. Therefore, this step 
should not aff'ect extended emission that has not already been 
filtered-out by the correlated noise removal. 

This reduction is optimized for recovering compact sources 
with high signal-to-noise ratio. Alternatively, the data can be re- 
duced without subtracting all components of correlated noise, 
which allows one to recover more extended emission, but results 
in an increased noise level in the final maps. 

3.2. Second step: using a source model 

The data processed in the first step are then combined into 
groups of three maps: for each position /q, the maps centered 
on /o, /o + 0.5°, and Iq - 0.5° are computed and coadded to obtain 
a two-fold coverage of the central 1x2 deg^ area (see Fig. O. 
This combined map is then used as a source model for mask- 
ing in the timestreams the data that correspond to positions in 
the map where the emission is above some given threshold. The 
full reduction is then performed again, involving the same pro- 
cessing steps as described in the previous section with less con- 
servative parameters for despiking and correlated-noise removal. 
During this step, datapoints corresponding to strong emission in 
the model map are masked for the various computations (base- 
line removal, despiking, subtraction of correlated noise), but the 
results of these computations are also applied to the masked data. 

Finally, another iteration is performed, this time subtracting 
the source model (resulting from the previous step) from the raw 
data. In this case, all datapoints are used in the computations, and 
processing steps that may aff'ect strong sources (e.g., despiking) 
can be performed almost without caution, since most of the true 
emission has been subtracted from the signals. Before adding 
the source model back into the reduced data, weights are com- 
puted using sliding windows along the timestreams. The weight 



applied to each datapoint is 1/rms^, where rms is the standard 
deviation in a given pixel computed for 50 datapoints. 

The net effect of these iterations with a source model is to 
reduce the negative artifacts appearing around strong sources in 
the first reduction step and to recover some fraction of the flux of 
bright objects lost in that step. The final flux uncertainty for com- 
pact sources should be lower than 15% (Sect. 13. iL as also con- 
firmed by h undreds of rneasure ments of primary and secondary 
calibrators (ISiringo et al.ll2009l) . 

When compiling maps with two pixels per beam, the fi- 
nal noise level obtained from individual observations is gener- 
ally in the 70-100 mJy/beam range. In the combined maps, the 
rms is then 50-70 mJy/beam in overlapping regions. The com- 
bined maps covering continuously -30° to -hi 2° in longitude are 
shown in Fig.[3j and example zooms are shown in Fig.lH 

3.3. Compact source extraction 

A preliminary extraction of co mpact sources was comp leted us- 
ing the SExtractor programme ( Berlin & Arnouts 1996). The ex- 
traction was performed in two steps: first, the peaks were de- 
tected on signal-to-noise ratio maps; SExtractor was then used 
to compute the source peak positions and fluxes, as well as their 
sizes. Since we were only interested in compact sources when 
compiling this first catalog, one constraint used was that the ra- 
tio of major to minor axes should be lower than 4. Details of 
the method and the complete catalog itself will be published in 
a forthcoming paper (Contreras et al. in prep). 

This source extraction resulted in -6000 objects with peak 
fluxes above 250 mJy/beam, over the 95 deg^ mapped. The aver- 
age source density is thus 63 sources per square degree, but the 
source distribution is highly non-uniform (Fig. [5]). The source 
sizes range from one LABOCA beam size (W) to about 8 beam 
widths (150''). We note that the filtering of uniform emission on 
spatial scales larger than 2f5 (Sect. 13.1b should not prevent ex- 
traction of sources of larger FWHM, as long as they are reason- 
ably compact, i.e., with a profile close to a Gaussian shape; their 
fluxes however would be severely underestimated. 

Assuming a typical distance of 5 kpc, which roughly cor- 
responds to the molecular ring, the 19" resolution translates to 
-0.5 pc. The LABOCA data are there fore sensitive to molec ular 
clu mps, following the terminol ogy of IWilliams et al.l (l2000l see 
also lMotte & Hennebelldl20 09V where dense cores have typical 
sizes of ~0.1 pc and are embedded in pc-scale clumps, them- 
selves being part of >10 pc GMCs. The sources distribution is 
strongly peaked toward the Galactic Plane, as seen in Fig.[5l and 
in particular in the direction of the Galactic Center, with an ex- 
cess of sources toward negative latitudes. Within the central few 
degrees, sources exhibit an excess at positive longitudes. Other 
peaks in longitude are visible in Fig. [5] in the direction of the 
Norma arm (-30° < / < -20°), and around / = -8°, which cor- 
responds to the NGC6334/NGC6357 star-forming complexes. 

As seen in Fig. [51 the latitude distribution of compact sources 
peaks at /? = -0.09°, which is significantly below zero. This ef- 
fect is seen in all ranges of longitudes. We suspect that this is due 
to the Sun being slightly above the Galactic Plane, but this eff'ect 
requires more study. The reason for the asymmetry in longitude 
within the few central degrees is unclear, but is certainly real, as 
can be seen also in Fig.[3l Interestingly, the distribution of com- 
pact sources seen with Spitzer a t 24 / im shows an excess toward 
negative longitudes (Hinz et al. 2009). While no clear explana- 
tion is found, Hinz et al. ( 2009) suggested that this could reflect 
the true distribution of infrared sources near the Galactic Center, 
if more molecular clouds are present at positive longitudes, ob- 
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Fig. 3. Combined maps showing 84 deg^ of the ATLASGAL data obtained in 2007, at 870 jim with a resolution of 19'\ The color 
scale is logarithmic between 0.05 and LO Jy/beam. The rms is of order 50 to 70 mJy/beam in most parts of these maps. A few 
known regions are labeled in the maps. See also Fig.lH 
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Fig. 4. Zooms into selected regions. These maps are projected with three pixels per beam, resulting in a pixel r.m.s. of between 80 
and 130 mJy/beam. All maps are shown in linear scale, calibrated to Jy/beam. The color scale has been chosen to highlight the 
negative artifacts, but doesn't cover the full range of fluxes in these images: the brightest peaks are, from left to right, 18.6, 17.2, 
and 7.2 Jy/beam; the most negative pixels are at -0.49, -0.54, and -0.43 Jy/beam, respectively. 
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Fig. 5. Distribution of the 6000 compact sources in Galactic lon- 
gitude (top, with 1° bins) and latitude (bottom, with 0.1° bins). 
Note that the longitude range -hll.5° to -hl5° is not covered by 
the present data. 



scuring the compact IR sources, which would be consistent with 
the distribution that we observe in the submm. 



4. First results 

The large scale maps in Fig. [3] clearly show a variety of submil- 
limeter sources from bright, compact objects to faint, extended 
regions on the arcmin scale, and long filaments on almost a de- 
gree scale. The survey is very rich, and a wealth of new results 
is expected. To provide a flavor of these and to illustrate how 
a complete view of the Galactic star-formation activity will be 
achieved from ATLASGAL, we present a first analysis and some 
first results for a typical 1° wide slice, around I = 19° (Fig. [6]). 
This region shows a reasonably rich population of sources, and 
is also representative of regions with several cloud complexes 
at diff'erent distances. We also selected this regi on becaiise it i s 
covered by the Galactic Ring Survey (GRS. Jackson et al.ll2006l) . 
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Fig. 6. Example combined map covering 1x1.5 deg^, shown in 
logarithmic scale. The pixel scale is one third of the beam and the 
standard deviation between pixels in empty regions of the maps 
is typically 90 mJy/beam. The brightest peak is at 6.2 Jy/beam. 
The 128 compact sources extracted in this region are shown with 
orange ellipses, with sizes and position angles as measured with 
SExtractor (Sect. 13.31) . Ultra-compact Hii regions discussed at 
the end of Sect. 14.21 are shown with blue ellipses; compact and 
classical Hii regions are indicated with magenta ellipses. 



which provides good quality CO data that we use to discuss the 
association with large-scale molecular clouds. 
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4.1. Compact sources and extended emission at I = 19° 

In a number of Galactic directions, there is submillimeter emis- 
sion virtually everywhere in the Galactic Plane. The region 
shown in Fig. [6] at / = 19° is one of these. From our prelimi- 
nary source extraction, we identified 128 compact sources with 
peak fluxes above 0.25 Jy/beam in this 1x2 deg^ slice (no source 
was detected outside the |/?| < 0.75° range shown in Fig. O, 
in perfect agreement with the average 63 extracted sources per 
square degree. Only 16 of these sources have a probable known 
counterpart in the SIMBAD database (infrared, millimeter or ra- 
dio source within a 10'' search radius), which means that more 
than 80% have not yet been studied. Some of them certainly 
have counterparts in the new infrared surveys observed with 
Spitzer (see e.g.. Fig. [14]); associations between ATLASGAL 
and Spitzer data will be addressed elsewhere. 

As can be seen in Fig. [51 most of the compact sources are 
embedded in more extended emission of a few arcmin extent. 
Many of these features appear filamentary and could correspond 
to overdensities within giant molecular clouds, in which even 
denser clumps are forming. That small groups of clumps be- 
long to one single molecular complex is sometimes confirmed 
by information about their kinematics, as discussed in the next 
section. 



4.2. Coherent complexes and distance determination 

Toward the inner Galactic Plane, the detected sources can be at 
any distance. A priori, any direction of the Plane could contain 
sources at distances distributed over a significant range of values, 
except for some peculiar lines of sight, like the Galactic Center 
(GC), where most sources are probably concentrated close to the 
GC itself. In other directions, one expects to find overdensities 
of objects at particular distances, corresponding to the Galactic 
arms, particularly in the molecular ring region, but a wide range 
of distances is allowed. 

Distance determination in the Galaxy is a difficult issue. 
The kinematic distances can be obtained from molecular-line 
follow-up surveys. We have started systematic observations of 
the compact ATLASGAL sources in NH3 with the Eff'elsberg 
100 m telescope (Wienen et al. in prep.) for the northern part of 
the plane, and in many molecular lines observed simultaneously 
with the MOPRA 22 m telescope (Wyrowski et al. in prep.) for 
the southern part. The NH3 lines were detected for nearly all tar- 
geted sources, showing various line widths and shapes, and usu- 
ally with only one velocity component (see examples in Fig. [7]). 
This is expected since NH3 is a high-density tracer, which is de- 
tectable only for the dense cores also seen with LAB OCA. 

Twenty-four of the total 278 sources observed so far in NH3 
are located in the region shown in Fig. [6] The distribution of 
measured Vlsr is shown in Fig. [S] Interestingly, the velocities 
are not spread over the entire permitted range, but are mostly 
gathered in very few coherent groups. In addition, these groups 
are also spatially coherent (see Fig.[8]top). We can conclude that 
the ATLASGAL emission is dominated by a few bright com- 
plexes at specific distances. In the / = 19° region, 15 out of 24 
ATLASGAL compact sources detected in NH3 have a velocity of 
between 60 and 66 km s"^ This complex would be at a distance 
of 4 or 12 kpc depending on whether the near or far solution for 
the kinematic distance is correct. In addition to this clear veloc- 
ity component, two sources are at rather low velocities of around 
20 km s"^ a few are between 35 and 50 km s"\ and two are at 
80 and 81 km s'^ 



To study the spatial distribution of these velocity components 
in greater detail, we used the ^^CO(l-O) data from the GRS to 
identify the molecular complexes associated with these veloci- 
ties (Fig. [9]). In addition, using an extinction map obtained from 
the 2MASS catalog (Bontemps, priv. com.), it is possible to rec- 
ognize that some of these complexes account for the detected 
extinction and are therefore most probably at their near kine- 
matic distance. This is the case for the 36 to 52 km s"^ complex, 
and for a small part of the -20 km s"^ clouds (dashed-line arc 
around (19.35, 0.05) in Fig. [9]). The main complex at 65 km s"^ 
is more complicated: parts of this complex are associated with 
peaks of extinction, but other parts correspond to very low ex- 
tinction. Nevertheless, this might indicate that it is also at the 
near distance. In addition, this complex appears quite extended 
at 8 yum and 870 yum, which again favors the near distance solu- 
tion. The full extent of the ^^CO emission, spanning about 0.9° 
in latitude, corresponds to 70 pc at a distance of 4.5 kpc, typical 
of GMC dimensions. 

Finally, we note that a large part of the ~20 km s"^ com- 
ponent is a string of compact CO clouds that are not associated 
with any feature in the extinction map. They are therefore most 
probably complexes at their far distance, i.e., on the order 14 
kpc. The source at (/,/?) = (18.7 6,0.26) coincides with a very 
compact object in the GLIMPSE (IB enj amin et al .1 [2QQ3b image 
(Fig. [To]), which is also consistent with a far distance solution. 

Additional constraint s can be obtained from absorption lines 
towar d Hii regions fe.g.JSewilo et aP l2004l: I Anderson & Banial 
120091) . lAnderson & Banial (l2009l) " presented seven Hii regions 
with measured distances determination in the range 18.5° < 
/ < 19.5°: three ultra-compact and four compact/classical re- 
gions. All three ultra-compact and three out of the four com- 
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Fig. 7. Example spectra in NH3 (1,1) transition for three sources 
located in the region shown in Fig.[6l G18.76-h0.26 is probably 
at a distance of -14 kpc, while G19.01-0.03 and G18.61-0.07 
are probably at 4.3 and 3.7 kpc, respectively. 
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Fig. 8. Top-. LSR velocities (in km s"^) derived from NH3 spectra 
are indicated next to the sources visible in the ATLASGAL map. 
Bottom', the measured LSR velocities are plotted against Galactic 
longitude. Kinematic distance s derived using the Galactic rota- 
tion model of Bran d & Blitzl 01993) are indicated on the right 
for the near distance solution, and on the left for the far solution. 
The thick line near 140 km s~^ shows the maximal expected ve- 
locity, which corresponds to the tangential point for a circular 
orbit with a radius ~2.7 kpc around the GC. 



pact regions have clear counterparts in ATLASGAL (Fig. [S]). 
Except for CI 8.95-0.02, which exhibits several components in 
the ^^CO spectrum, and C19. 04-0.43, which we did not observe 
in NH3, there is good agreement between the LSR velocities de- 
rived from radio recombination lines, from |^CQ(l-0) and from 
NH 3 transitions. Finally, resul ts from, both JSewilo et al.l (|2004|) 
and I Anderson & Banial (l2009l) generally agree with our conclu- 
sions for the objects located in the test field, the -h65 km s"^ 
complex being at the near distance and the -h20 km s"^ clump 
G19.47-h0.17 being at the far distance. However, in a few 
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Fig. 9. Integrated ^^CO(l-O) intensity over diff'erent ranges of 
velocities, shown in the top-left corner of each panel, from the 
GRS data. The color scale is shown on top of each panel, in 
units of K km s"^. LSR velocities of ATLASGAL sources de- 
rived from the NH3 spectra are indicated. Bottom right: extinc- 
tion map derived from the 2MASS catalog. The color scale cor- 
responds to magnitudes of visual extinction. In all panels, dashed 
lines outline structures that seem connected together. 



cases, the correct interpret ation of absorption lines remains dif- 
ficult: ISewilo et al.l (l2004 and I Anderson & Banial (l2009h con- 
cluded with a far distance solution for Ul 8.66-0.06, but all three 
ATLASGAL sources that we observed in NH3 (with vlsr ~ 
-h45 km s"^) follow an IRDC clearly detected by GLIMPSE 
(Figs. [TOl and [14]), which favors a near distance. Discussion of 
all individual sources goes beyond the scope of the present pa- 
per, and will be addressed elsewhere. 

4.3. Mass estimates 

Dust emission is generally optically thin in t he submm contin- 
uum. Therefore, following Hildebrandl (Il983|) . the total mass in 
a core is directly proportional to the total flux density Fy, inte- 
grated over the source: 



M = 



d^FyR 

Bv(Td) Ky 



(1) 



where d is the distance to the source, R is the gas-to-dust mass 
ratio. By is the Planck function for a dust temperature To, and 
Ky is the dust absorption coefficient. The beam-averaged column 
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Fig. 11. Total gas plus dust masses in a core with 250 mJy inte- 
grated flux density at 870 jim (5-cr detection limit) as a function 
of dust temperature, for three typical distances: 500 pc (full line), 
2 kpc (dashed line) and 8 kpc (dotted line). 
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Fig. 10. 8 yum map of the region shown in Fig. [6] from the 
GLIMPSE survey. LSR velocities derived from NH3 data are 
indicated. 

density, which does not depend on the distance, can also be com- 
puted with: 



By{TD)^KyiimH 



(2) 



where Q is the beam solid angle, ji is the mean molecular weight 
of the interstellar medium, which we assume to be equal to 2.8, 
and m// is the mass of an hydrogen atom. Assuming a gas-to- 
dust mass ratio of 100, and Ky = L85 cm^ (interpolated to 
870 yum from Table 1, Col. 5 of Ossenkopf & Henning|Ll994), 
our survey is sensitive to cold cores with typical masses in the 
range -100 M© at 8 kpc to below 1 M© at 500 pc (Fig. fTTT) . 
The 5-cr limit of 0.25 Jy/beam corresponds to column densities 
in the range 3.6 x 10^^ cm'^ for Td = 30K to 2.0 x 10^^ cm'^ for 
T/) = 10 K, equivalent to visual extinctions of 4 mag to 21 mag, 
respe ctively (using the conversion factor from iFrerking et al.l 
fT982h . 

These computations assume that the emission at 870 yum is 
optically thin. For the brightest source detected so far 
(Sgr B2(N), which has a peak flux density at 170 Jy/beam), we 
compute the beam-averaged optical depth using: 



T870 = -In 



1 - 



^ByiT) 



(3) 



This infers an optical depth of r < 0.31 for a temperature T > 
25 K. Therefore, the optically thin assumption should be correct, 
even for the brightest sources. 

Using the distances as estimated in Sect. 14.21 we calculated 
the masses of the clumps located around / = -hl9° (Fig. [8]); the 



results are summarized in Table [T] In a few cases, the SExtractor 
program, used to extract compact sources from the full survey 
(Sect. [33]), does not reliably extract sources that appear blended. 
Therefore, we fitted 2D-Gaussian profiles to the 24 sources ob- 
served in NH3 to measure their fluxes and sizes, special care 
being taken to separate blended sources. The results of these fits 
were source positions, sizes (minor and major axis FWHM), and 
both peak and integrated fluxes. 

Table [T] provides: source names (based on Galactic co- 
ordinates); equatorial coordinates and sizes, measured from 
Gaussian fits on the 870 yum image; LSR velocities measured 
with the NH3(1,1) transition; corresponding Galactocentric ra- 
dius, as well as near and far kinematic distances; peak and in- 
tegrated fluxes at 870 yum; and kinetic temperatures, computed 
from the NH3(1,1) and (2,2) measurements (Wienen et al., in 
prep.). Finally, masses for the near and far distances are com- 
puted, using the integrated 870 yum fluxes and the kinetic tem- 
peratures. 

These computations assume that the 870 yum fluxes detected 
with LAB OCA originate only in thermal emission from the 
dust. From spectral-line surveys, it is known that molecular lines 
might contribute a significant am ount of flux to the emission 
detected by broadband bolometers. iGroesbeckl (1 19951) estimated 
from line surveys in the 345 GHz atmospheric window that lines 
can contribute up to 60% of the flux. This extreme value was 
found for the Orion hot core, w hich is peculiar due to a large con- 
tribution from wide SO2 lines fechil ke et al.lll997h . The relative 
contribution from CO(3-2) lines alone will only be of relevance 
in the case of extreme outflow sources with little dust contin- 
uum. From a quantitative analysis of the relative CO contribu- 
tion in the bright massive star-forming r egion G327. 3-0.6, ob - 
served with the APEX telescope in CO ( Wvrowsk i et al.ll2006h . 
we found the contribution to be insignificant close to the com- 
pact hot core but to increase to 20% for the bright, extended 
photon-dominated region to the north of the hot core. In sum- 
mary, hot molecular cores, strong outflow sources, and bright 
photon-dominated regions will have a flux contribution from line 
emission that is higher than the typical 15% calibration uncer- 
tainty (Sect. 13. ll) in the continuum observations. Line follow-up 
observations in the mm/submm with APEX of bright sources 
found with ATLASGAL will help to place stronger constraints 
on the line flux contribution. 

Another source of contamination could be free-free emis- 
sion, when an Hii region has already been formed. At 345 GHz 
the contribution from free-free should however be low compared 
t o the thermal emission from dust. As an example, Motte et alj 
(i2003h estimated the contribution from free-free emission at 
1.3 mm in the compact fragments found within the W43 massive 
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Table 1. Observed and derived parameters for the 24 sources located in the region shown in Fig.[6]and observed in NH3 lines. When 
the distance ambiguity could be solved, the preferred distance and mass are shown in boldface. 



Source 


RA (J2000) 


Dec (J2000) 


Size 


Vlsr 


Rg 


Dnear 


Dfar 


-Ppeak 


Pint 


Tkin 


N(H2) 


-^^near 


Mfar 










km/s 


kpc 


kpc 


kpc 


Jy/b. 


Jy 


K 


(1) 


Mo 


Mo 


G18.55+0.04 


18 24 37.85 


-12 45 14.8 


36x24 


36.3 


5.6 


3.1 


13.0 


0.92 


2.8 


21.0 


2.1 


140 


2380 


G18.61-0.07 


18 25 08.42 


-12 45 21.6 


38x26 


45.8 


5.2 


3.7 


12.5 


1.70 


5.7 


32.7^2) 


2.2 


210 


2490 


G18.63-0.07 


18 25 09.84 


-12 44 07.5 


53x21 


45.7 


5.2 


3.6 


12.5 


0.96 


3.6 


14.7 


3.9 


420 


5010 


G18.65-0.06 


18 25 10.63 


-12 42 24.1 


43x25 


44.6 


5.2 


3.7 


12.5 


2.25 


8.1 


18.3 


6.4 


680 


7940 


G18.66+0.04 


18 24 50.98 


-12 39 20.8 


36x24 


80.9 


4.0 


5.2 


10.9 


2.10 


6.3 


22.4 


4.5 


780 


3510 


G18.70+0.00 


18 25 02.38 


-12 38 12.6 


46x34 


80.3 


4.0 


5.1 


11.0 


0.84 


4.5 


14.8 


3.4 


1030 


4770 


G18.76+0.26 


18 24 13.27 


-12 27 44.6 


59x42 


20.8 


6.6 


2.0 


14.1 


1.14 


9.7 


15.8 


4.1 


320 


15090 


G18.79-0.29 


18 26 15.53 


-12 41 36.2 


40x34 


65.5 


4.4 


4.5 


11.6 


0.94 


4.3 


14.4 


4.0 


830 


5360 


G18.83-0.30 


18 26 23.70 


-12 39 38.9 


37x28 


63.5 


4.6 


4.4 


11.7 


1.35 


4.8 


_(3) 


5.3 


800 


5610 


G18.83-0.47 


18 26 59.34 


-12 44 46.5 


34x32 


63.3 


4.5 


4.5 


11.6 


1.50 


5.7 


18.9 


4.1 


680 


4600 


G18.83-0.49 


18 27 03.57 


-12 45 06.5 


58x25 


65.2 


4.5 


4.5 


11.6 


1.34 


6.7 


20.8 


3.2 


720 


4640 


G18.86-0.48 


18 27 06.20 


-12 43 09.3 


38x23 


65.4 


4.4 


4.6 


11.5 


0.73 


2.1 


18.8 


2.0 


270 


1710 


G18.88-0.49 


18 27 09.21 


-12 42 35.1 


38x25 


65.5 


4.5 


4.5 


11.6 


1.11 


3.5 


22.4 


2.4 


340 


2210 


G18.88-0.51 


18 27 15.13 


-12 42 54.6 


52x29 


66.2 


4.4 


4.6 


11.5 


0.62 


3.2 


23.6 


1.2 


290 


1840 


G18.89-0.47 


18 27 07.82 


-12 41 38.0 


64x34 


66.1 


4.4 


4.6 


11.5 


3.39 


25.4 


25.9 


5.9 


2040 


12840 


G18.89-0.64 


18 27 43.46 


-12 45 48.5 


28x22 


64.1 


4.5 


4.5 


11.6 


0.58 


1.2 


20.0 


1.5 


130 


880 


G18.91-0.63 


18 27 43.55 


-12 44 51.0 


50x43 


64.5 


4.5 


4.5 


11.6 


1.51 


11.1 


17.1 


4.8 


1560 


10340 


G18.93-0.01 


18 25 32.35 


-12 26 41.4 


27x25 


46.3 


5.1 


3.7 


12.4 


1.33 


3.1 


29.3 


2.0 


140 


1530 


G18.98-0.07 


18 25 51.06 


-12 25 42.5 


26x18 


61.8 


4.6 


4.4 


11.7 


0.82 


1.3 


20.4 


2.0 


140 


960 


G18.98-0.09 


18 25 55.22 


-12 25 50.3 


30x19 


63.0 


4.5 


4.4 


11.6 


0.64 


1.3 


20.7 


1.5 


130 


900 


Gl 8.99-0.04 


18 25 45.44 


-12 23 54.0 


38x32 


59.8 


4.7 


4.3 


11.8 


0.31 


1.3 


20.7 


0.7 


120 


930 


G19.01-0.03 


18 25 44.61 


-12 22 42.0 


40x34 


59.8 


4.7 


4.3 


11.8 


2.18 


10.1 


19.5 


5.7 


1070 


8000 


G19.08-0.29 


18 26 48.77 


-12 26 20.7 


39x37 


66.0 


4.6 


4.3 


11.7 


4.28 


21.4 


24.2 


8.2 


1700 


12330 


G19.47+0.17 


18 25 54.50 


-11 52 34.4 


33x24 


20.0 


6.6 


2.0 


14.0 


6.53 


18.0 


23.4 


13.1 


330 


15510 


Notes: Unit of N(H2) is 10^^ 


-2 

cm 

























This source was also observed in NH3 bv lSridharan et al.l j2005h . with a better signal-to-noise ratio. We therefore use the kinetic temperature 



derived by this study rather than our own data for this source. 

Only the (1,1) transition of NH3 is detected at this position, so that no kinetic temperature could be derived. A temperature of 15 K has been 
assumed to compute masses and column density for this source. 




Flux density [Jy] 

Fig. 12. Distribution of peak fluxes for the compact sources ex- 
tracted in the data observed in 2007, shown in logarithmic scale. 
Fractions of sources associated with IRAS or MSX sources 
within a 3C search radius are indicated, using the scale shown 
in the right hand side. 



star-forming region, which harbors a giant Hii region. In most 
cases, the free-free emission accounts for <40% at 1.3 mm, one 
extreme case being where this fraction reaches 70%. Assuming 
a flat spectrum for free-free (optically thin case), and Fy oc y^-^ 
for the dust emission (i.e., a dust emissivity of j3 = 1.5), scal- 
ing from 240 GHz to 345 GHz implies that 20% of the flux 
originates from free-free emission at the LABOCA wavelength. 
In conclusion, the contribution from free-free should be almost 
negligible, and even in the most extreme cases should not exceed 
20%. 

4.4. Nature of the sources 

Since dust emission in the submm is optically thin, the 
LABOCA observations are sensitive to all kinds of objects con- 



taining dust, which can range from cold, dense clumps within 
clouds to more evolved proto-stars that are centrally heated. We 
searched systematically for infrared counterparts to our prelimi- 
n ary catalog of 6000 s ources in the IRAS Point Source Catalog 
(Beichman et al. and in the Midcourse Space Experiment 

(MSX, Price et al. 2001) data, using a search radius of 30'' 
(more details will be given in Contreras et al. in prep.). Only 
-10% of the submm compact sources can be associated with an 
IRAS point source, and -32% have a possible MSX counter- 
part (Fig.[T2l). In contrast, two thirds do not have a mid- or far- 
infrared counterpart in the IRAS/MSX catalogs, and could cor- 
respond to the coldest stage of star formation, prior to the birth of 
a massive proto-star. This fraction is similar to the results found 
by Motte et al. (2007) for the Cygnus-X region. However, with 
the superior sensiti vity of the Spitzer GLIMPSE and MIPSGAL 
(ICarey et al.l l2005l) surveys, weaker mid-IR sources are found 
to be associated with many of the MSX-dark, compact submm- 
clumps. Sources with IRAS associations generally correspond to 
later stages of star formation, e.g., hot molecular cores, compact 
Hii regions, or young embedded star clusters. 

One should note that in some bright, complex star-forming 
regions, because of crowding and blending, even strong infrared 
sources do not fulfill the point- source selection criteria and, con- 
sequently, are not included in the IRAS/MSX catalogs. An ex- 
treme example is seen in the histogram of Fig. [121 where the 
brightest source in the survey, SgrB2(N), is counted as a source 
without IRAS/MSX association. 

High-contrast IRDCs are easily identified by ATLASGAL. 
These clouds can have various shapes, from compact cores to 
filamentary morphologies of 1-5 arcmin in size. As an example. 
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Fig. 13. Column density map of the IRDC G19.30+0.07, con- 
tained in the I = 19° demonstration field. The map is derived 
from 8 Jim extinction values measured in the corresponding 
Spitzer/GLIMPSE map. The ATLASGAL emission at 870 yum is 
overlaid in contours. The map center is at Of = 18^ 25"^ 53.94^, 
S = -12° 05^43.9'' (J2000), or / = 19.2111% b = +0.0703°. 



our demonstration fi eld contains the already well-investigated 
IRDC G19.30+0.07 dCarev et al.l Il998). In Fig. [13 we show 
a column-density map of this IRDC, based on an extinction 
map derived fr om GLIMPSE 8 jum d a ta. Th e Ry = 5.5 ex- 
tinction law of IWeingartner & Draini (1200 ll) has been used. 
Further details of this approach are described in a forthcom- 
ing publication (IVasyunina et al .1120091) . Overlaid as contours are 
the ATLASGAL 870 yum emission data. There is good agree- 
ment between the morphologies observed in both datasets, espe- 
cially at the high column-density center of the IRDC. With an 
adopted distance of 2.38 kpc, roughly 260 Mq are collected in 
this clump. IRDCs are considered to be a good hunting ground 
for the earliest phases of star formation, and ATLASGAL is ex- 
pected to reveal a considerable fraction of them to be sub-mm 
emission sources. 

Some sources seen in the submm may also correspond to 
low-mass star-forming sites located at small distances. However, 
our survey is sensitive to 1 Mq dense clumps up to a distance of 
~1 kpc only (Eq. (1)). In addition, nearby regions of star for- 
mation are likely to be located away from the Galactic Plane, 
outside the area covered by ATLASGAL. Therefore, low-mass 
star-forming clumps are not expected to be a dominant popula- 
tion in our data. 

A small fraction of sources detected in the submm may not 
be related to star formation. These sources include late evolved 
stars with circumstellar envelopes, planetary nebulae, the nuclei 
of external galaxies, and quasars. However, we estimate that the 
contamination due to these classes of obj ects is negligib l e. For 
instance, following the method outlined bylAnglada e t al.l (1 19981 
Appendix), the number of background sources per arcmin^ with 
flux densities at 5 GHz above S is: 

NiarcmirT^) = 0.01 1 X ^ (mJyY^-^^ (4) 

Assuming that the flux density scales with frequency as 5" oc 
y"^-^, one then obtains, for y = 345 GHz, 0.10 sources above 
150 mJy per square degree, or about 10 background sources in 
the 95 deg^ area discussed here. The extrapolation from 5 GHz 
to 345 GHz is quite uncertain, but this example illustrates nev- 



ertheless that the number of these background sources must be 
very low. 

Another estimate can be derived from deep-field sur- 
veys in the submm. Integrating the diff'erential source 
counts from the S CUBA Half-Degree Extragalactic Survey at 
850 yum rSHADES. ICoppin et al.l2006l their Eq. 12 and Table 7), 
yields 1 x 10"^ sources per square degree with flux densities 
above 150 mJy, or 0.1 sources above this limit in our surveyed 
area. Thus, we do not expect to detect any high-redshift object 
in our data. 

The situation for late evolved stars may be more complex. 
As an example, the red supergiant VY CMa, which is one of 
the brightest sources at IR wavelengths and is at a distance of 
-1.5 kpc, has a flux density of 1.5 Jy at 870 yum, as measured 
with LAB OCA. Therefore, this object is detectable out to a dis- 
tance of 3.5 kpc with our survey sensitivity. Less extreme or 
more distant objects of similar nature will most likely fall below 
our detection limits, but we cannot exclude having a small frac- 
tion of supergiants and extreme late evolved stars in our data. It 
will be possible to identify such objects from their colors at other 
wavelengths, e.g., in the Spitzer surveys (Robitaille et al.„20Q8h . 

5. Perspectives 

5.1. An unbiased survey 

This survey will provide the first unbiased database of Galactic 
sources at submm wavelengths, covering the inner 120° of the 
Galaxy in a systematic fashion. Since dust emission in the 
submm is optically thin, it is directly proportional to the mass of 
material (assuming a constant temperature) and leads to the de- 
tection of all proto-stellar objects or pre-stellar clumps and cores. 
With a 5-(T detection limit of 250 mJy/beam, this survey will be 
complete for proto-stellar cores down to ~10 Mq at 2 kpc, or 
-50 M© at 5 kpc (Fig.[TT]). In addition, the data will allow the 
investigation of the structure of molecular complexes on vari- 
ous scales, from their lower-intensity filaments to higher-density 
star-forming regions. 

The systematic nature of this survey will also allow a com- 
plete census of extreme sources in the inner Galaxy, such as the 
high-mass star-forming regions W49 or W5 1 . With our sensitiv- 
ity, we should detect, at a 10-cr level, sources on the far side of 
the Galaxy, 25 kpc away, of masses as low as about 1500 Mq. 

The unbiased nature of this database is critically important 
in compiling robust samples of objects for all types of follow-up 
projects based on a reliable statistical analysis. This is especially 
true in the context of ALMA (see also Sect. 15. 51 below), which 
will explore exactly the same part of the sky as APEX in the 
(sub)mm range. 

5.2. Evolutionary stages of high-mass star-formation 

This database will contain a large number of high-mass clouds 
with submm emission not associated with developed star for- 
mation (as traced in the infrared surveys), which are expected 
to correspond to the earliest phases during the formation of the 
richest clusters of the Galaxy. Since high-mass proto-stars are 
rare (iMotte et al.ll2007l) . such large-scale surveys are required to 
obtain large samples of objects corresponding to the diff'erent 
stages of massive-star formation. Thus, the data will allow us 
to determine a basic evolutionary sequence for high-mass star 
formation. 

Another topic that can be addressed with these data is the 
importance of triggered star formation, which can be assessed 
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by identifyi ng the second generatio n stars at the edges of Hii re- 
gions (e.g., Deharveng et al. '2005) and supernova remnants. In 
the context of high-mass star formation in particular, the impor- 
tance of triggering compared to spontaneous star formation re- 
mains an open question. 

5.3. Distances and Galactic structure 

Dust emission relates the star-forming regions seen in the 
IR/radio to the molecular clouds seen in CO, e.g., in the Galactic 
Ring Survey, by tracing the highest column-density features: 
compact sources and filaments. With the help of both extinc- 
tion m aps and absorption lines toward Hii regions (ISewilo et al.l 
|2004|) . the near/far ambiguity of kinematic distances can be 
solved in deriving the proper distances for a large fraction of 
ATLASGAL sources (Sect. l4~2l) . This will therefore help to im- 
prove our 3D view of the Galaxy. These data will also allow us to 
discover new coherent star-formation complexes up to distances 
of a few kpc. Finally, these data will allow us to compare the 
star-formation processes at diff'erent locations, in the spiral arms 
and inter-arm regions, and between the diff'erent arms. 

5.4. Ancillary data 

The region that we plan to map covers the same area as the 
GLIMPSE and MIPSGAL surveys (e.g., Eig. [H). It will also 
be mapped with PACS and SPIRE for the approved Hi-GAL 
(the Herschel infrared Galactic Plane Survey; Molinari et al.) 
open-time Herschel key project. The ATLASGAL data will 
thus be cross-correlated with similar large-scale Galactic sur- 
veys at other wavelengths, in the near- to far-infrared range, 
for example with the /, H, K UKIDSS Galactic Plane Survey 
(iLucas et al.ll200 8), providing spectral energy distributions from 
~1 to 870 /im. The LABOCA data at 870 fim are crucial to es- 
timating the total amount of dust along the line of sight. The 
combination of data from APEX, Spitzer, and Herschel will pro- 
vide an extraordinary new database of the Galactic Plane from 
the near-IR to the millimeter, at a spatial resolution ^30''. 

We will also correlate our data with existing radio surveys 
in the centinieter ra nge: the NRAO VLA Sky Survey (NVSS, 
ICondon et al.lll998|). the Mult i-Arrav Galactic Plane Imaging 
Survey (MAGPIS, |Helfeniet al. 2006) , the Sydn ey University 
Molonglo Sky Survey (SUMSS, Mauch et al.ll2003l) . and the Co- 
ordinated Radio 'N' Infrared Survey for High-mass star forma- 
tion (CORNISnfl It will be of particular interest to compare 
with surveys of the radio continuum and maser emission that 
will be possible with the Expanded Very Large Array and th e 
Australia Telescope Compact Array in the near future (iMentenI 
120071) . Existing or upcoming molecular-line databases include 
those resulting from the multi-beam CH3OH maser Galactic 
Plane survey (Green et al. 2007) a nd from surveys, e.g., in the 
CS o r ^^CO molecules (Bronfman et al.| Il996| : [Urquhart etaP 
I2OO7I) . Combining ATLASGAL with the IR to cm data will allow 
us to establish a solid picture of the star-formation evolutionary 
stages and their associated timescales. 

Several bolometer arrays are expected to be installed at 
APEX in the near future, including the Submillimetre APEX 
Bolometer Camera (SABOCA), a 37-element array observing 
at 350 yum, with a beam of 7fl. Our team plans to complement 
the 870 yum data with targeted observations at 350 yum, cover- 
ing limited areas of a few arcmin extent, to constrain the dust 
emissivity and resolve bright sources in crowded regions. The 
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Fig. 14. Three-color image from Spitzer surveys: GLIMPSE 
4.5 yum in blue (mostly tracing stars) and 8 yum in green (mostly 
PAH emission), MIPSGAL 24 yum in red, with 870 yum emission 
overlaid in contours. 



870 fim survey might also be extended in surface area (toward 
the outer Galactic Plane), or to greater depth for a small part of 
the Galactic Plane (20 mJy/beam rms or less). 

A similar survey of the northern Galactic Plane (/ > -10°) 
has started with the BOLOCAM instrument at the Caltech 
Submillimeter Observ atory (CSO), at a wavelength of 1.1 mm 
(iDrosback et al.l l2008h . with a 27'' resolution. The northern 
Galactic Plane will also be mapped with the SCUBA-2 instru- 
ment at 450 and 850 yum, with resolutions of 8'' and 15'', re- 
spectively. In particular, as part of the SCUBA-2 "All-Sky" 
Survey (SASSY. Thompson et al.ll2007h . a 10° wide strip cov- 
ering \b\ < 5° over the part of the Galactic Plane accessible 
from Hawaii will be observed at 850 yum with a one-cr sensitiv- 
ity of 30 mJy/beam, comparable to the ATLASGAL survey. In 
addition, the JCMT Galactic Plane Survey plans to cover much 
deeper the area -hlO° < / < -h65°, \b\ < 1°, at both 450 and 
850 yum. Combining data from these large-scale surveys will pro- 
vide data at three diff'erent wavelengths in the Rayleigh- Jeans 
part of the dust emission spectrum, thus allowing us to study, 
among others, the dust emissivity, related to the composition of 
the dust. 

5.5. Legacy value 

The raw data will be available to outside users through the ESO 
archiveEl one year after the end of each observing period. We 
also plan to make the reduced combined maps available to the 
community; in addition, a compact- source catalog (including 
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associations with surveys at other wavelengths) and a catalog 
of extended objects will be published upon completion. These 
will be available from a public web siteEI. This project will thus 
have a high legacy value, providing well characterized samples 
of targets to address in detail many questions of astrophysics. 
Calibrated maps built from the data observed in 2007 and dis- 
cussed in the present paper are expected to be available online, 
both on the ESO archive and the ATLASGAL public website, by 
Summer 2009. 

The observations that we have planned should be completed 
by the end of 2009, i.e. when the Herschel observatory is ex- 
pected to deliver its first science data. Our team will be in- 
volved in the investigation of systematic associations between 
the ATLASGAL sources and Herschel products. In particular, 
members of the consortium are leading the key projects HOB YS 
(the Herschel imaging survey of OB Young Stellar objects; 
Motte, Zavagno, Bontemps et al.) and Hi-GAL. 

The first ALMA antennas should also be available for early 
science soon after completion of our survey. This will be the per- 
fect time for conducting follow-ups of carefully selected targets 
at an unprecedented spatial resolution. The ATLASGAL cata- 
log will provide ALMA with more than 10^ compact sources, 
including hot cores, high-mass proto-stellar objects, and cold 
pre- stellar cores. All sorts of molecular lines will be available 
to study the chemistry of the sources, or to search for signa- 
tures of infall or outflows for example. Thus, ATLASGAL is 
expected to trigger numerous follow-up studies to characterize 
the detected sources in terms of their excitation and chemistry. 
ATLASGAL will also provide the James Webb Space Telescope 
(JWST) with samples of well characterized targets, and these fu- 
ture facilities will allow one to study the inner workings of the 
star-formation process on small scales. This survey is therefore 
an essential pathfinder for Herschel, ALMA, and the JWST. 

In addition to the new findings in the realm of high-mass star 
formation, the survey will tell us on a larger scale details about 
the structure and mass distribution in our Galaxy and help us to 
connect our knowledge of local star formation with star forma- 
tion in external galaxies, where high-mass stars are the only ones 
that we see. 
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